ANALYSIS  OF  THE  HF  DATA  MEASURED  AT 
LEITRIM,  USING  A  MATRIX  CALIBRATION 
METHOD  DEVELOPMENT  AT  DREO  (U) 


by 

Eric  K.L.  Hung  and  Zhijiang  Hu 


fcwo  WAUrr  arasBOM d  a 


DEFENCE  RESEARCH  ESTABLISHMENT  OTTAWA 


Canada 


REPORT  NO.  1319 

"DiaTRiSlJTION  STArLteMT  A 

Approved  for  public  release; 
Distribution  Urriiinit§d 


19971224  042 


November  1997 
Ottawa 


1+1 


National  Defense 

Defence  nationale 


ANALYSIS  OF  THE  HF  DATA  MEASURED  AT 
LEITRIM,  USING  A  MATRIX  CALIBRATION 
METHOD  DEVELOPMENT  AT  DREO  (U) 


by 

Eric  K.L.  Hung 

Ground  Based  Group 
Surface  Radar  Section 

and 

Zhijiang  Hu 

Telexis  Corporation 


DEFENCE  RESEARCH  ESTABLISHMENT  OTTAWA 

REPORT  NO.  1319 


PROJECT 

05B09 


November  1997 
Ottawa 


ABSTRACT 


This  is  the  second  report  on  the  analysis  of  the  calibration  data 
taken  with  a  high  frequency  (HF)  antenna  array  and  a  signal  source  towed 
behind  a  small  aircraft.  The  analysis  described  in  this  report  uses  the 
second  of  two  array  calibration  methods  developed  by  E.K.L.  Hung  at  DREO; 
the  previous  report  uses  the  first  method.  This  work  is  part  of  a  bigger 
task  to  develop  a  calibration  method  for  the  HF  antenna  arrays  in  the 
Canadian  Forces  Supplementary  Radio  Systems  (CFSRS)  modernization 
project . 

The  calibration  method  used  in  this  report  is  superior  to  the  method 
used  in  the  previous  report.  The  results  derived  from  the  calibration 
data  generated  with  source  elevations  31°  or  lower  show  that  it  reduces 
the  deviation  of  the  array  beam  patterns  from  the  ideal  ones,  reduces  the 
errors  in  signal  direction  estimates,  and  increases  the  ability  of  a 
high-resolution  direction-finding  method  to  resolve  signal  directions. 


RESUME 


Le  present  rapport  est  le  second  sur  1' analyse  des  donnees 
d'etalonnage  recueillies  a  l'aide  d'un  reseau  d'antennes  a  haute 
frequence  (HF)  ainsi  que  d'une  source  de  signaux  trainee  par  un  petit 
avion.  Dans  ce  rapport,  1' analyse  est  realisee  au  moyen  de  la  seconde  des 
deux  methodes  d'etalonnage  du  reseau,  mis  au  point  par  E.K.L.  Hung  du 
CRDO ;  le  rapport  precedent,  quant  a  lui,  traitait  de  la  premiere  de  ces 
methodes.  Cette  analyse  s'inscrit  dans  une  tache  plus  vaste  portant  sur 
la  mise  au  point  d'une  methode  d'etalonnage  pour  les  reseaux  d'antennes 
a  haute  frequence  dans  le  cadre  du  pro jet  de  modernisation  du  Reseau 
radio  supplementaire  des  Forces  canadiennes  (RRSFC) . 

Dans  ce  rapport,  la  methode  d'etalonnage  est  superieure  a  celle  que 
decrit  le  rapport  precedent.  Les  resultats  decoulant  des  donnees  de 
calibrage,  qui  proviennent  de  sources  ayant  une  elevation  d'au  plus  31°, 
indique  une  reduction  de  la  deviation  entre  les  diagrammes  de  faisceau 
du  reseau  et  les  donnees  ideales,  une  reduction  des  erreurs  dans  les 
evaluations  de  direction  de  signal,  et  une  augmentation  de  la  capacite 
de  la  methode  radiogoniometrique  a  haute  resolution  qui  sert  a  trouver 
la  direction  des  signaux. 
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EXECUTIVE  SUMMARY 


This  report  contains  a  second  analysis  of  the  calibration  data  taken 
with  a  high  frequency  (HF)  antenna  array  in  1994  at  Leitrim,  near  the 
Ottawa  airport.  The  analysis  uses  the  second  antenna  array  calibration 
method  developed  by  Eric  K.L.  Hung  at  DREO.  (The  first  report  is  on  the 
analysis  of  the  first  method.)  It  is  part  of  a  bigger  task  to  develop  a 
calibration  method  for  the  HF  antenna  arrays  in  the  Canadian  Forces 
Supplementary  Radio  Systems  (CFSRS)  modernization  proj ect .  These  arrays 
are  used  for  detecting  long-range  HF  signals  and  estimating  their 
directions . 

A  perfectly  calibrated  antenna  array,  in  conjunction  with  a  high- 
resolution  direction- finding  (HRDF)  method,  will  have  the  lowest  signal- 
to-noise  ratio  (SNR)  requirement  for  target  detection  and  the  most 
accurate  direction  estimates  for  the  array  and  HRDF  method.  Calibration 
errors  will  raise  the  minimum  SNR  required  for  target  detection  and 
reduce  the  accuracy  of  the  target  direction  estimates. 

The  method  here  calibrates  an  antenna  array  through  the  construction 
of  a  calibration  matrix  for  the  antenna.  It  differs  from  the  other  matrix 
methods  published  in  the  open  literature  in  that  it  is  tolerant  to  the 
presence  of  interfering  signals. 

The  following  observations  are  based  on  the  analysis  of  data  files 
generated  with  signal  source  elevations  31°  or  lower: 

1.  The  calibration  method  reduces  the  deviation  of  the  antenna  beam 
patterns  from  the  ideal  patterns; 

2.  The  method  increases  the  accuracy  of  direction  estimates  and 
suppresses  the  dependence  of  errors  in  bearing  estimates  on  source 
bearings;  and 

3  .  The  method  increases  the  ability  of  a  HRDF  method  to  resolve  signal 
directions . 

4.  The  method  is  superior  to  the  first  antenna  array  calibration 
method.  In  particular,  the  direction  estimates  are  more  accurate  and 
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the  ability  of  an  HRDF  method  to  resolve  signal  directions  is 
better . 

It  is  not  possible  to  remark  on  the  results  obtained  with 
calibration  data  generated  with  signal  source  elevations  higher  than  31°. 
The  results  are  very  encouraging.  However,  these  data  files  have  large 
measurement  errors  that  cannot  be  suppressed. 
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1 . 0  INTRODUCTION 


The  evaluation  of  an  antenna  array  calibration  method  [1]  developed 
at  DREO  is  presented  here.  It  uses  the  calibration  data  taken  in  1994 
with  a  high  frequency  (HF)  antenna  array  at  Leitrim,  which  is  near  the 
Ottawa  International  Airport.  The  array  received  signals  from  a 
transmitter  towed  behind  a  small  aircraft.  The  objective  was  to  develop 
a  calibration  method  for  the  HF  antenna  arrays  in  the  CFSRS  (Canadian 
Forces  Supplementary  Radio  System)  Modernization  Project. 

High-resolution  direction-finding  (HRDF)  methods  are  used  to  process 
the  outputs  of  HF  antenna  arrays  in  the  detection  of  targets  and  the 
estimation  of  their  directions.  To  fully  exploit  the  capabilities  of 
these  methods,  the  arrays  must  be  accurately  calibrated.  Only  in-situ 
calibration  methods  can  be  used,  because  the  arrays  are  physically  very 
large.  Ideally,  the  methods  should  account  for  the  mutual  coupling  among 
the  array  elements,  the  interaction  between  the  elements  and  local 
objects,  small  errors  in  element  positions,  and  the  presence  of 
interfering  signals. 

An  effective  way  to  calibrate  an  antenna  array  is  to  move  a  signal 
source  over  a  fine  grid  of  directions  covered  by  the  array  and  measure 
the  array  response  vectors  (array  manifold  vectors)  for  these  directions. 
However,  this  method  is  impractical,  because  the  computational  load  in 
direction  finding  with  these  vectors  is  very  high.  Besides,  different 
sets  of  array  response  vectors  are  needed  for  direction  finding  in 
differtent  environmental  conditions  (such  as  wet  spring  and  dry  autumn) . 

The  alternative  methods  treat  antenna  array  calibration  as  a 
parameter  estimation  problem.  They  note  that  the  outputs  of  the  array 
elements  are  usually  not  ideal  and  are  corrupted  by  error  sources  such 
as  unequal  element  gain  and  phase  responses,  mutual  coupling  among  the 
array  elements,  and  small  errors  in  element  positions.  These  methods 
calculate  a  set  of  calibration  coefficients  for  the  array  and  use  them 
to  suppress  the  errors  at  the  outputs  of  the  array  elements. 

Many  parametric  methods  for  calibrating  narrowband  antenna  arrays 
have  been  published  in  the  open  literature.  They  include  those  by  Paulraj 
and  Kailath  [2]  ,  Friedlander  and  Weiss  [3,4]  ,  Brown,  McClellan  and  Holder 
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[5,6],  and  Pierre  and  Kaveh  [7],  Hung  [1,8],  See  [9,10],  and  Dogan  and 
Mendel  [11] .  All  methods  account  for  unequal  element  gains  and  phases. 
The  method  in  [8]  differs  from  the  others  in  that  it  is  the  only  method 
that  can  also  account  for  the  presence  of  specular  multipath 
interference.  This  method  has  been  evaluated  by  Huang  and  Hung  in  [12] , 
using  the  above  calibration  data  collected  at  Leitrim.  The  methods  in 
[1]  /  [4]  ,  [7]  and  [9]  and  [10]  account  for  mutual  coupling  and 
approximate  its  effect  by  a  matrix  operation. 

The  matrix  method  in  [1]  is  evaluated  here,  because  it  accounts  for 
mutual  coupling  and  preliminary  studies  indicated  that  it  is  also  highly 
tolerant  to  the  presence  of  interfering  signals.  The  other  methods  that 
also  account  for  mutual  coupling  are  ignored  here.  There  are  several 
reasons.  Firstly,  Hung  has  shown  [13]  that  the  method  in  [4]  sometimes 
has  non-unique  solutions.  Secondly,  the  method  in  [7]  produces 
calibration  matrices  that  are  highly  dependent  on  the  calibration  data. 
The  reasons  are  given  by  Hung  [1]  and  See  [9] .  Thirdly,  the  methods  in 
[9]  and  [10]  assume  the  absence  of  interfering  signals.  This  assumption 
is  usually  invalid  in  the  calibration  of  HF  arrays. 

The  results  of  this  evaluation  show  that  the  method  in  [1]  reduces 
the  deviation  of  array  beam  patterns  from  the  ideal  beam  patterns, 
reduces  the  errors  in  direction  estimates,  and  increases  the  ability  of 
an  HRDF  method  to  resolve  signal  sources  close  in  direction. 

In  report  [12]  ,  it  was  also  noted  that  the  method  in  [8]  reduces  the 
errors  in  direction  estimates  and  increases  the  ability  of  an  HRDF  method 
to  resolve  signal  directions.  When  the  figures  produced  by  these  methods 
are  compared,  we  note  that  the  method  in  [1]  is  superior.  In  particular, 
the  errors  in  direction  estimates  are  smaller  and  do  not  depend  on  source 
bearings  .,  Besides,  with  the  method  in  [1]  ,  the  ability  of  an  HRDF  method 
to  resolve  signal  directions  is  slightly  better. 

The  remainder  of  this  report  is  organized  as  follows:  Chapter  2 
contains  a  description  of  the  antenna  array  and  the  coordinate  systems 
used  in  the  analysis.  Chapter  3  describes  the  calibration  measurements 
and  the  data  collected.  Chapter  4  discusses  the  preprocessing  of  these 
calibration  data;  it  includes  a  check  on  the  possible  presence  of  faulty 
array  elements,  Doppler  filtering,  correction  for  the  deviation  of  signal 
wavefronts  from  spherical  wavefronts,  and  the  reduction  of  the  sizes  of 
data  files.  The  matrix  calibration  method  is  presented  in  Chapter  5,  and 
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the  calibration  matrices  it  produces  are  discussed  in  Chapter  6.  In 
Chapter  7  is  an  evaluation  of  the  calibration  method.  It  studies  the 
antenna  array  beam  patterns,  the  accuracy  of  direction  estimates,  and  the 
ability  of  a  HFDF  method  to  resolve  signal  directions.  The  conclusions 
are  given  in  Chapter  8 . 
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2 . 0  ANTENNA  ARRAY 


The  antenna  array  was  located  at  Leitrim,  which  is  near  the  Ottawa 
International  Airport.  It  was  planar,  had  twelve  elements,  and  was  part 
of  a  larger  array  with  52  elements.  In  Table  1  and  Figure  1  are  the 
element  positions.  Given  below  are  the  element  parameters: 

Elements  1  to  8  - 

2  inches  diameter  aluminum  pipe,  vertical,  20  feet  high, 
antenna  feed  7  feet  above  ground,  eight  8-feet  ground  radials. 

Elements  9  to  12  - 

2  inches  diameter  aluminum  pipe,  vertical,  20  feet  high, 
antenna  feed  7  feet  above  ground,  no  ground  radial. 

Elements  9  to  12  were  installed  in  the  middle  of  winter,  years  after 
elements  1  to  8 . 

There  are  two  coordinate  systems  in  Table  1:  a  UVW- rectangular 
system  supplied  with  the  data  and  a  polar  system  to  be  used  in  the 
analysis  of  calibration  data.  These  systems  are  defined  in  Figure  2.  In 
the  rectangular  system,  the  U-direction  is  the  north  direction,  the  V- 
direction  is  west,  and  the  W-direction  is  up.  In  the  polar  system,  the 
location  of  a  point  is  identified  by  a  triplet  (r ,d,\p).  Here,  r  is  the 
distance  of  the  point  from  the  origin,  6  is  the  bearing  angle,  and  \f/  is 
the  elevation  angle.  The  bearing  6  is  measured  in  the  UV-plane  and 
clockwise  from  the  U-axis. 


The  rectangular  coordinates  (u,v,w)  and  polar  coordinates  (r ,d,\p) 
of  a  point  in  space  are  related  to  each  other  by: 


and 


u  =  r  cos \p  cosd  , 
v  =  -  r  cos\p  sinfl  , 
w  =  r  sini/'  , 

r  =  (u  +v  +w  )  , 

6  =  tan’1  ( -v/u)  , 

\p  =  tan’1  (w/  (u2+v2) 1/2) 


(2.1) 

(2.2) 

(2.3) 

(2.4) 

(2.5) 

(2.6) 
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The  bearing  6  is  related  to  the  azimuth  angle  6 az  in  an  XYZ-coordinate 
system,  defined  with  {X, Y, z} = {East , North, Up} ,  by 

0az  =  90°  -  6  .  (2.7) 

The  coordinates  in  the  UVW  and  XYZ  systems  are  related  by 


x  =  -v  , 
y  =  u  , 
z  =  w  . 


(2.8) 

(2.9) 

(2.10) 


N 

Rectal 

u 

ngle  Coord 

V 

ELEMENT  P 

inate 

w 

OSITIONS 

Pole 

r 

ir  Coordin 

8 

ate 

1 

3 .263 

-24 .786 

0 . 0 

25 

82 . 5 

0 . 0 

2 

-3 .263 

-24 . 786 

o 

o 

25 

97 . 5 

o 

o 

3 

-15 .219 

-19 . 834 

o 

o 

25 

127 . 5 

o 

o 

4 

-19 . 834 

-15 .219 

o 

o 

25 

142 . 5 

0 . 0 

5 

-24 .786 

3 .263 

0 . 0 

25 

187.5 

0.0 

6 

-19 . 834 

15.219 

0 . 0 

25 

217 . 5 

0 . 0 

7 

23 . 097 

9.567 

0 . 0 

25 

337 . 5 

0 . 0 

8 

24 . 786 

3.263 

0 . 0 

25 

352 . 5 

0 . 0 

9 

85.001 

22 . 776 

o 

o 

88 

345 . 0 

o 

o 

10 

-22 . 776 

85 . 001 

o 

o 

88 

255 . 0 

o 

o 

11 

120 . 741 

32.352 

o 

o 

125 

345 . 0 

o 

o 

12 

-32 .352 

120 . 741 

0 . 0 

125 

255 . 0 

0 . 0 

Table  1.  Element  positions  in  the  HF  antenna  array.  Distances 
u,  v,  w,  and  r  are  measured  in  meters.  Angles  8  and  i p  are 
measured  in  degrees . 
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3 . 0  CALIBRATION  MEASUREMENTS 


The  calibration  data  were  measured  with  a  vertically  oriented  dipole 
transmitter  towed  behind  an  airplane.  This  plane  was  flown  in  specified 
paths  while  snapshots  of  the  outputs  of  the  array  elements  were  taken. 
Two  global  positioning  receiver  systems  were  used  to  determine  the  ranges 
and  directions  of  the  transmitter  in  the  snapshots.  The  first,  identified 
as  GPS1 ,  was  at  the  antenna  site.  Its  output  time,  denoted  by  tl,  was 
recorded  with  the  array  snapshots.  The  second,  identified  as  GPS2,  was 
a  differential  system  consisting  of  two  GPS  receivers,  one  on  the  plane 
and  the  other  at  a  known  ground  location.  Its  output  time  t2  and  position 
rGPS2  ( t2 )  were  recorded  in  a  log  file  on  the  aircraft  at  approximately  one 
second  intervals. 

Six  transmitter  frequencies  were  used:  5.1,  7.5,  9.3,  11.5,  15.1, 
and  18.0  MHz.  The  outputs  at  each  array  element  were  downshifted  in 
frequency  to  near-baseband  with  commonly  tuned  receivers,  sampled  with 
an  A/D  converter  at  10.4k  samples/second,  and  stored  in  blocks  of  8000 
array  snapshots.  Approximately  fifty  data  files  were  generated. 

This  study  is  restricted  to  the  analysis  of  the  data  files  generated 
with  a  transmitter  frequency  of  5.1  MHz  only.  In  these  files,  the  near¬ 
basebands  are  centred  at  an  intermediate  frequency  (IF)  of  approximately 
2.0  kHz  and  the  filter  bandwidths  at  the  mixers  were  3.0  kHz. 

Table  2  contains  the  nominal  flight  parameters  and  the  sizes  of  the 
files  at  5.1  MHz.  Figures  3  to  10  illustrate  the  paths  of  the  airplane 
in  these  files.  They  include  five  clockwise  (CW)  circles  round  the  array, 
one  counter-clockwise  (CCW)  circle  round  the  array,  and  two  straight 
lines  over  the  array. 

The  flight  path  in  File  2  has  a  side  loop  near  bearing  210°  and  a 
hook  near  bearing  333°.  In  this  study,  we  refer  to  File  2  as  the  modified 
file  with  the  side  loop  and  the  hook  removed. 

The  transmitter  was  towed  behind  the  airplane  and  transmitted  a 
continuous  wave  signal  to  the  array.  Its  location  at  the  recorded  GPS2 
time  t2  was  approximately  the  same  as  that  of  GPS2  at  time  t2-A,  where 
A  was  equal  to  the  GPS2-to-transmitter  distance  divided  by  the  speed  of 
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the  airplane.  The  value  of  A  was  approximately  2  seconds;  its  exact  value 
depended  on  the  wind  speed,  wind  direction,  plane  speed,  and  plane 
direction  at  time  t2 . 

In  a  separate  study  of  the  calibration  data  by  Jenkins  [14]  ,  it  was 
noted  that  the  GPS1  time  tl  and  the  GPS2  time  t2,  recorded  at  the  same 
time  instant,  were  related  by 

t2  =  tl- €  ,  (3.1) 

where  e  was  a  non- zero  difference  that  was  essentially  independent  of 
frequency  and  transmitter  elevation.  The  origin  of  this  non- zero 
difference  was  unknown.  The  study  also  showed  that  the  transmitter 
position  at  time  tl,  denoted  by  rTX(tl) ,  could  be  derived  from  the  GPS2 
position  at  t2  as 

rTx  ( ^  1 )  =  rGPS2  ( 1 2  -  A ) 

—  ^gps2  (tl-A),  (3.2) 

where  A'=e+A.  The  assignment 

A'  =  8.7  seconds  (3.3) 

produced  the  best  overall  agreement  between  the  estimated  and  measured 
transmitter  directions. 
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File 

Date 

Elevation 

(deg) 

Plannee 

Altitude 

(ft) 

Track 

Radius 

(km) 

Direction 

Data 

Blocks 

Recorded 

1 

18/05/94 

4 . 5 

3000 

10 

CW 

2255 

2 

20/05/94 

11 

5000 

7 

CW 

2308 

3 

20/05/94 

11 

5000 

7 

ccw 

1438 

4 

24/05/94 

18 

6000 

5 

CW 

1106 

5 

24/05/94 

31 

6150 

3 

CW 

723 

6 

24/05/94 

49 

6000 

2 

CW 

367 

7 

24/05/94 

6000 

1350 

8 

24/05/94 

6000 

1365 

Table  2.  Nominal  flight  parameters  for  the  data  files  measured  at  5 . 1 
MHz.  The  date  format  is  day/month/year .  Each  data  block  contains  8000 
array  snapshots. 
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Figure  3.  Flight  path  in  File  1.  The  airplane  circles  the 
array  with  a  radius  of  10  km  and  proceeds  clockwise  (CW) 
through  405°,  from  north  45°  east  to  north  90°  east. 

U(km) 


Figure  4 .  Flight  path  in  File  2  before  the  removal  of  the 
side  loop  at  bearing  210°  and  hook  at  bearing  at  333°.  The 
airplane  circles  the  array  with  a  radius  7  km,  proceeds  CW 
through  358°,  from  north  24°  west  to  north  26°  west. 


U(km) 


Figure  5.  Flight  path  in  File  3:  range  7  km;  counter-CW  (CCW) 
405°,  north  19°  west  to  north  64°  west. 


U(km) 


Figure  6.  Flight  path  in  File  4:  range  5  km;  CW  404°,  south 
49°  west  to  north  57°  west. 
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Figure  9.  Flight  path  in  File  7.  It  is  a  straight  line  over 
the  array  with  a  heading  of  west  9.7°  north. 


U(km) 


Figure  10.  Flight  path  in  File  8.  It  is  a  straight  line  over 
the  array  with  a  heading  of  north  1.0°  east. 


4.0  DATA  PREPROCESSING 


The  data  files  require  preprocessing  before  they  can  be  used  to 
evaluate  the  calibration  method  in  [1]  .  This  preprocessing  comprises  the 
following  tasks: 

1.  Check  if  there  were  faulty  elements  during  the  experiment; 

2.  Doppler  filter  the  data  in  each  file  to  increase  the  signal- to-noise 
ratios  (SNRs)  in  the  array  snapshots  by  30  dB  and  to  reduce  the 
number  of  snapshots  in  the  file  by  a  factor  of  1000; 

3.  Reduce  the  size  of  each  filtered  file  further  to  approximately  3600 
snapshots  with  approximately  0.1°  spacing  in  signal  source  bearings; 
and 

4.  Convert  the  signal  wavefronts  in  the  reduced  files  to  planar 
wavefronts . 

Task  2  is  necessary,  because  the  array  snapshots  in  the  original 
files  have  low  SNRs  and  there  are  2.9  to  19  million  snapshots  per  file. 
Task  4  is  also  necessary,  because  the  transmitter  ranges  (track  radii) 
are  too  short  in  the  files. 

There  are  two  measurement  errors  that  cannot  be  corrected  in  the 
preprocessing.  The  first  is  intrinsic  to  the  antenna  array  and  originated 
from  the  misalignment  of  the  array  elements.  It  has  been  assumed  that  the 
aluminum  pipes  used  in  the  array  elements  were  vertical  and,  therefore, 
parallel.  The  relative  orientations  of  some  pipes  were  checked  with  the 
parallax  method  in  the  summer  of  1995.  They  were  not  exactly  parallel, 
as  the  angle  between  two  pipes  was  several  degrees  in  some  cases.  This 
misalignment  perturbed  the  gain  and  phase  responses  of  the  array 
elements.  The  perturbation  increased  with  transmitter  elevations  and 
varied  slowly  with  the  transmitter  bearing.  The  second  error  is  an 
artifact  that  originated  from  the  short  transmitter  ranges  in  data 
generation.  These  ranges  were  10  km  or  less.  Consequently,  one  cannot 
ignore  the  difference  in  the  transmitter  elevations  at  the  array 
elements.  This  error  also  perturbed  the  gain  and  phase  responses  of  the 
elements.  Usually,  the  perturbation  was  larger  if  the  transmitter  range 
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was  shorter. 


File  6  in  Table  2  was  measured  with  a  rather  high  transmitter 
elevation  of  49°  and  a  rather  short  transmitter  range  of  2.0  km.  The 
measurement  errors  due  to  the  above  perturbations  are  so  large  that  this 
file  is  not  expected  to  produce  reliable  results  in  the  study.  Files  7 
and  8  have  too  many  snapshots  measured  with  transmitter  ranges  shorter 
than  2  km  and  elevations  higher  than  49°.  They  are  not  used  in  this 
study. 

4 . 1  Faulty  Element  Check 

By  design,  the  transmitter  signal  at  the  array  elements  was  very 
strong  compared  with  the  other  signals.  The  frequency  of  this  signal  in 
the  data  files  is  close  to  the  IF  (=2.0  kHz)  at  the  receivers,  because 
the  nominal  transmitter  range  was  constant  and  the  corresponding  Doppler 
frequency  of  the  transmitter  was  zero.  Therefore,  the  antenna  array  was 
unlikely  to  have  faulty  elements  if  (a)  the  spectrum  of  the  outputs  of 
each  element  has  a  very  strong  signal  present,  (b)  the  frequency  of  this 
signal  is  the  same  for  all  elements,  and  (c)  this  frequency  is  close  to 
2 . 0  kHz . 

Some  results  obtained  in  a  check  for  faulty  elements  are  presented 
in  Figs.  11  and  12.  They  were  produced  with  the  Fast  Fourier  Transform 
(FFT)  and  the  first  1000  array  snapshots  in  the  first  block  of  Files  1 
and  2,  respectively.  Each  figure  has  twelve  spectra,  one  from  the  outputs 
of  each  array  element.  Because  the  figures  are  crowded  and  the  spectra 
in  the  same  figure  look  alike,  the  element  identification  numbers  have 
been  omitted.  A  separate  study  shows  that,  in  each  spectrum,  the 
strongest  amplitude  has  a  frequency  of  1.95  kHz  and  is  at  least  25  dB 
stronger  than  the  other  signals.  Therefore,  the  array  was  unlikely  to 
have  faulty  elements  during  the  generation  of  Files  1  and  2. 

Three  identifiable  interfering  signals  are  present  in  the  spectra. 
Their  frequencies  are  0.00,  0.06  and  3.08  kHz.  The  0.00  Hz  signal 
probably  came  from  offsets  in  the  A/D  converters.  The  origins  of  the 
others  are  unknown. 

Similar  studies  showed  that  the  antenna  array  was  also  unlikely  to 
have  faulty  elements  during  the  generation  of  Files  3  to  6 . 
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Figure  12 .  The  spectra  of  the  outputs  of  the  twelve  array  elements  in 
Block  1  of  File  2 . 
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Only  100  0  array  snapshots  were  used  to  calculate  each  Doppler 
spectrum,  although  each  data  block  had  8000  snapshots.  This  number  of 
snapshots  was  selected,  because  the  bin  width  in  the  FFT  spectrum  must 
not  be  smaller  than  the  frequency  bandwidth  of  the  received  transmitter 
signal  and  the  Doppler  shift  due  to  transmitter  motion  was  changing  with 
time . 


To  produce  a  spectrum,  the  sequence  of  1000  real  output  amplitudes 
of  an  array  element  was  augmented  with  1048  zeros;  the  augmented  sequence 
was  converted  to  a  2048  point  power  spectrum,  using  the  FFT;  and  the 
first  1024  amplitudes  in  the  spectrum  were  plotted.  The  augmentation  with 
zeros  was  necessary.  Otherwise,  picket  fence  effect  and  noise  presence 
could  combine  to  prevent  the  strongest  amplitude  from  appearing  at  the 
same  frequency  bin  in  every  spectrum.  Only  the  first  1024  amplitudes  were 
plotted,  because  the  element  output  amplitudes  were  real  and  did  not  have 
quadrature  components.  Consequently,  the  spectrum  was  symmetrical  about 
bin  k=l ,  and  the  amplitudes  for  bins  k  and  (2050-k)  were  the  same. 

Each  FFT  spectrum  has  2048  frequency  bins.  These  bins  are  identified 
with  k=l  to  2048.  The  frequency  at  the  k-th  bin  is  5.0781(k-l)  Hz  for  k=l 
to  1024  and  5 . 0781 (k-2049)  Hz  for  k=1025  to  2048.  It  is  equal  to  | f s- f 0 1 
where  fs  is  the  signal  frequency  and  fo  is  the  oscillator  frequency. 

4.2  Doppler  Filtering 

Doppler  filtering  was  used  to  reduce  each  data  block  from  8000  array 
snapshots  to  eight  new  ones  3  0  dB  stronger  in  SNR  compared  with  the 
original  snapshots.  The  procedure  is  given  below. 

Step  1. 


Separate  the  data  in  each  data  block  into  eight  sets  with  1000 
consecutive  snapshots  per  set. 

Step  2  . 

For  each  set, 

(a)  Augment  the  1000  output  amplitudes  with  1048  zeros; 

(b)  Calculate  the  2048  point  complex  FFT  spectrum  of  the  augmented 
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sequence;  and 


(c)  Retain  the  first  1024  complex  amplitudes. 


Step  3  . 

Study  the  complex  FFT  spectrum  calculated  with  the  outputs  of  the 
first  array  element.  Identify  the  signal  bin  as  the  bin  with  the 
largest  power. 


Step  4. 


Construct  a  new  array  snapshot  as 

vll  _  /  -y  II  -y  II  y  1 1  ...  y  II  ^ 

A  —  \  -A-  x  2  '  ^  3  ;  ^  12' 

Denote  the  corresponding  transmitter  location  by  r"  . 


(4.1) 


In  (4.1)  ,  x1}  to  x'^2  are  the  complex  amplitudes  in  the  signal  bins  in 
the  FFT  spectra  calculated  with  the  outputs  of  elements  1  to  12, 
respectively.  The  superscript  'T'  denotes  the  transpose  and  the  location 
r",  where  r"=(r "  ,6"  ,\p")  ,  denotes  the  average  of  the  transmitter  locations 
in  the  1000  array  snapshots. 


4.3  File  Size  Reduction 


This  task  is  designed  to  reduce  the  size  of  each  filtered  file  to 
M  snapshots,  with  M=3600  for  Files  1  to  5  and  M=2435  for  File  6.  File  6 
has  a  smaller  M,  because  it  has  less  than  3000  snapshots  before 
reduction.  In  the  following  procedure,  the  set  of  array  snapshots  before 
reduction  is  denoted  by  { ,  x"2 ,  *"3 ,  •  •  •  }  and  the  set  after  reduction  is 
denoted  by  { x'x ,  x'2 ,  x'3 ,  ■  •  •  ,  x'„}  . 


Step  1 . 

Generate  a  set  of  M  dummy  bearings  as 

6  =  0,  +  360 (m-1) /M  .  m=l , 2 , ■ ■ ■ , M-l  ,  (4.2) 

m  1 

where  0X  is  the  bearing  of  the  first  snapshot  in  the  filtered  file. 
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Step  2  . 


In  the  set  {x'J,x''  ,x'^,  •  •  •  }  ,  identify  the  snapshot  with  bearing 
closest  to  0m  as  x'm.  Denote  the  transmitter  location  of  this 
snapshot  by  r'm. 

The  r'm  in  Step  2  represents  the  set  (r'm,  6'm,  ,  where  r'm  is  the 

measured  transmitter  range. 

4.4  Correction  for  Non- Planar  Wavefront 

Antenna  arrays  are  normally  calibrated  to  receive  signals  from 
remote  sources.  These  sources  produce  planar  signal  wavefronts  at  the 
array.  However,  the  wavefronts  in  {x; , x'2 , x'3 ,  •  ■  •  , x'J  were  not  planar, 
because  the  transmitter  range  was  10  km  or  less.  Given  x'  (snapshot  x'm 
with  identification  suppressed)  ,  we  correct  the  deviation  from  a  planar 
wavefront  by  replacing  it  with  a  new  snapshot  x,  where  x=(xx,x2,  •  ■  •  ,x12)T, 
constructed  as 

xn  =  {(sn/rt)  exp[j27r(sn-rt)/A]  x'J  a n(0t,^t)  , 

n=l,2, • • • , 12.  (4.3) 


In  this  expression, 

{s1(  s2,  s3,  •  •  •  ,  s12)  =  element  distances  from  transmitter, 

rt  =  measured  transmitter  range  from  origin  of  coordinate 
system, 

A  =  signal  wavelength 
=  58.82  m, 


and 


=  exp  [ j 2tt  (uncosi/'tcos0t  -  vncosi/'tsinet 

+  wnsini/'t)  /A]  .  (4.4) 

In  (4.3),  the  factor  (sn/rt)  corrects  the  r'2  dependence  of  signal 
power  on  range,  exp  [  j27r  (sn-rt) /A]  equalizes  the  signal  phases  at  the 
elements,  and  an(8t,\pt)  is  the  n-th  component  of  the  array  steering  vector 
for  direction  {6t,\pt)  .  In  the  absence  of  noise  and  interference,  the 
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product  inside  the  bracket  {  }  is  the  same  for  all  array  elements  and  the 
source  direction  of  x  is  determined  by  an(6t,\pt)  ,  n=l  to  12.  The  signal 
wavefront  is  planar  in  x. 
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5.0  CALIBRATION  METHOD 


A  brief  description  of  the  calibration  method  in  [1]  is  presented 
here.  It  includes  the  signal  model  used,  the  assumptions  made,  and  the 
procedure  for  calculating  the  calibration  matrix. 

In  the  presentation,  N  denotes  the  number  of  array  elements,  the  set 
{l-g2-g3/ •  •  • -gN}  denotes  the  element  gains,  { 0 , 62,  <53 ,  •  •  •  ,  6N}  denotes  the 
element  phases,  {x1,x2,x3,-  •  •  ,xM}  denotes  M  preprocessed  array  snapshots, 
and  {  (61,\p1)  ,  (62,\p2)  ,  {d3,\p3)  ,  •  •  ■  ,  }  denotes  the  signal  directions  in 

the  snapshots.  The  coefficients  of  a  matrix  B,  for  example,  are  denoted 
by  (kPq:  p = 1 , 2  ,  ■  •  •  ;  q=l ,  2  ,  ■  •  •  }  .  B  is  a  normalized  matrix  if  bxl=l. 

The  element  gains  and  phases  at  array  elements  2  to  N  are  unknowns 
and  are  measured  relative  to  those  of  the  first  element.  They  are 
independent  of  direction.  Their  deviations  from  unit  gain  and  zero  phase 
originate  from  mismatches  in  the  cables  and  electronics  at  the  elements. 

5.1  Signal  Model 

In  the  signal  model  used,  the  set  of  M  preprocessed  array  snapshots 
are  associated  with  a  set  of  M  ideal  array  snapshots  {y3 ,  y2 ,  y3 ,  •  •  •  ,  yM) 
given  by 


yra  =  m=l ,  2  ,  •  •  •  ,M  ,  (5.1) 

where  am  is  the  complex  signal  amplitude  in  snapshot  ym,  and  vm  is  the 
array  steering  vector  for  the  corresponding  source  direction  .  This 

expression  assumes  a  perfectly  calibrated  array,  a  pure  tone  signal,  no 
interaction  among  the  array  elements,  no  noise,  and  no  interference. 

The  measured  array  snapshots  are  related  to  the  ideal  snapshots  by 

xm  =  TB(ym+em),  m=l,2,-**,M  ,  (5.2) 

where  r  is  the  matrix  of  element  gains  and  phases  and  is  given  by 

r  =  diag  { 1 ,  g2exp  [  j  62]  ,  g3exp  [  j  63]  ,  •  •  •  ,  gNexp  [  j  6J  }  ,  (5.3) 
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B  is  a  normalized  matrix  related  to  the  coupling  among  the  array 
elements,  and  em  is  an  error  term  representing  the  sum  effect  of  noise, 
interfering  signals,  sensor  misplacements,  interactions  between  the  array 
elements  and  local  objects,  etc.. 

The  calibration  matrix  is  defined  as  the  normalized  matrix  C  that, 
together  with  a  set  of  amplitude  estimates  {o^ ,  a2 ,  a3 ,  •  •  •  ,  aM}  ,  minimizes  the 
ob j  ective  function 

M 

S  =  I  |Cxm-yJ2  .  (5-4) 

m=l 


i 


.  e . 


{C,a}  =  arg  min  { S }  .  (5.5) 

C,a 

Function  S  depends  on  the  unknowns  to  aH  through  (5.1)  .  The 
denotes  the  estimate.  For  convenience,  a  has  been  used  to  denote  the  set 
{5'i,Qf2,Q'3,  '  •  '  ,«„}  • 

Because  C  is  normalized,  it  satisfies  the  condition 

clx  =  1.0  .  (5.6) 

This  condition  ensures  that  C  is  unique.  It  also  excludes  the  trivial 
solution  {c=0N,  q'1=0!2='  •  *  =o;H=0 }  ,  where  0N  denotes  the  NxN  null  matrix. 

5.2  Assumptions 

The  following  assumptions  are  made  in  the  development  of  the 
calibration  method: 

1.  The  array  elements  are  narrowband; 

2.  The  mutual  coupling  among  the  elements  can  be  approximated  by  a 
matrix  operation; 

3.  The  signal  source  (the  transmitter)  moves  over  a  wide  range  of 
directions  during  the  calibration  experiment  ,- 
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4.  The  number  of  array  snapshots  in  the  calibration  data  is  reasonably 
large ; 

5.  A  preliminary  calibration  has  been  carried  out  to  equalize  the 
element  gains  and  phases  of  the  antenna  array.  After  this 
calibration,  the  global  maximum  of  every  antenna  beam  pattern  is  the 
maximum  closest  to  the  beam  steering  direction;  and 

6.  The  signal-to-interference-plus-noise  ratio  (SINR)  in  the  array 
snapshots  is  high  enough  for  estimating  the  complex  signal 
amplitudes  in  the  snapshots. 

5.3  Computation  of  Calibration  Matrix 

The  calibration  matrix  C  is  obtained  by  iteration,  assuming  that 
only  one  signal  source  is  present  in  every  array  snapshot.  In  each 
iteration,  the  sum  S  is  initially  minimized  with  respect  to  the  estimates 
of  signal  amplitudes.  Then  it  is  minimized  with  respect  to  C.  The 
procedure  is  based  on  the  following  observations: 

A.  If  Ckl  is  an  estimate  of  C  and 


z 


tn 


(5.7) 


the  least  square  estimate  of  cvm  is  given  by 

~  H  /  I  |  2 

“m  =  VmZm/lVJ  ' 

and  the  corresponding  estimate  of  ym  is  given  by 


(5.8) 


y  =  O'  V 
«•  m  mm 


(5.9) 


B.  Given  {Yj , y2 ,  •  ■  ■  , yM}  ,  one  can  get  a  new  normalized  least  square 
estimate  of  C  by  solving  the  equation 


CkX  =  ,Y 


(5.10) 


as 


U  ^  tJ  _  1 

ck  =  (yy”)  (XY") 


(5.11) 
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and  then  impose  the  condition  (5.6)  on  Ck . 

In  these  expressions,  X  denotes  the  NxM  matrix  (xlfx2,  •  •  •  ,  xH)  ,  Y 
denotes  (y1 , y2 ,  •  •  •  , yM)  ,  and  the  superscript  'H'  denotes  the 
conjugate- transpose .  The  inverse  (XY11)"1  exists,  because  of 
Assumption  3 . 

C.  The  initialization  C0=IN,  where  IN  is  an  NxN  identify  matrix, 
is  equivalent  to  using  the  original  data  in  the  first 
iteration,  where  k=l. 


D.  Let  a  difference  matrix  be  defined  as 


Let  the  magnitude  of  the  largest  coefficient  in  this  matrix  be  given 
by  ek,  i  .  e  .  , 


ek  =  max {  |  dpq  |  :  p,  q=l ,  2  ,  •  •  •  ,  N} 


(5.13) 


Where  each  d  is  a  coefficient  of  Dk.  One  can  terminate  the 
iteration  if  (a)  ek<e' ,  where  e'  is  a  prespecified  small  value  or 
(b)  ek>ek.1 .  In  case  (a),  the  iteration  has  converged  to  a 
prespecified  level.  In  (b)  ,  the  convergence  condition  is  not 
satisfied  yet,  but  more  iterations  do  not  guarantee  more  accurate 
estimates  of  C. 


In  a  preliminary  study  with  File  2  and  100  iterations,  ek  decreased 
steadily  from  0.44064  for  k=l  to  0.00035  for  k=23 .  It  then  varied 
randomly  between  0.00016  to  0.00063  for  the  next  77  iterations.  With  File 
6,  ek  decreased  steadily  from  2.70607  for  k=l  to  0.06215  for  k=14.  It  then 
varied  randomly  between  0.02910  and  0.15553  in  the  next  86  iterations. 

The  computation  of  the  calibration  matrix  is  described  below.  There 
are  six  steps. 


Step  1 . 


Let  k=l  and  C0=IN. 
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Step  2 . 


Calculate  {y1,y2,  •  •  •  ,  yM }  with  (5.7)  to  (5.9). 

Step  3  . 

Calculate  Ck  with  (5.11) 

Step  4 . 


Impose  condition  (5.6)  on  Ck . 

Step  5 . 

Calculate  ek  with  (5.12)  and  (5.13). 


Step  6  . 

If  ek<0.0005  or  ek>ek.1, 

terminate  the  iterations  and  identify  Ck  as  C. 
Else, 

replace  k  by  k+1  and  return  to  Step  2. 


A  few  notes  are  in  order.  Firstly,  in  Step  2  of  the  last  iteration, 
to  aH  are  estimates  of  the  complex  signal  amplitudes  in  x1  to  xM, 
respectively.  Secondly,  given  a  calibration  matrix  C  and  a  measured  array 
snapshot  x,  the  corresponding  calibrated  array  snapshot,  denoted  by  xc, 
is  given  by 

Xc  =  Cx  .  (5.14) 

Thirdly,  from  (5.2)  and  (5.4),  the  normalized  inverse  of  C,  denoted  by 
C’1,  is  related  to  the  product  TB  by 

c'1  =  b^(rB)  ,  (5.15) 

where  blx  is  the  first  diagonal  element  of  B.  Finally,  if 


26 


bn  =  b22  =  '  '  ‘  =  bNN  (Assumption)  , 

one  can  estimate  r  and  the  product  b^B  as 

f  =  matrix  C  1  with  off-diagonal  elements 
equated  to  zero  , 


(5.16) 

(5.17) 

(5.18) 
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6.0  CALIBRATION  MATRICES 


Each  one  of  the  six  preprocessed  data  files  has  been  used  to 
calculate  a  calibration  matrix  C,  a  normalized  inverse  C1  of  C,  an 
estimate  of  r,  and  an  estimate  of  the  product  b^B.  Due  to  the  hugh 
quantity  of  data  generated,  only  the  notable  observations  are  presented 
here.  In  the  discussions,  C(2),  for  example,  denotes  the  matrix 
calculated  with  File  2.  The  corresponding  estimate  of  b‘*B  is  denoted  by 
bf*B(2)  . 


The  assumption  (5.16)  is  only  approximately  valid  here .  However,  the 
results  obtained  with  (5.17)  and  (5.18)  are  included  in  the  discussions, 
because  the  estimates  f1  and  b^B  contain  special  features  that  can  be 
identified  with  the  antenna  array  configuration  and  its  environment. 

The  preprocessed  data  files  are  used  in  the  calculation,  because 
they  have  many  desirable  properties .  Firstly,  in  a  preliminary  study,  the 
SINRs  in  the  individual  array  snapshots  are  usually  over  15  dB .  Secondly, 
the  number  of  snapshots  per  file  is  not  too  large  and  not  too  small. 
Thirdly,  the  signal  source  (the  transmitter)  circles  the  array  once  in 
the  file.  Fourthly,  the  snapshots  are  approximately  equally  spaced  in 
source  bearings . 


The  estimates  of  element  gains  and  phases  are  in  Tables  3  and  4, 
respectively.  They  were  calculated  with  (5.3)  and  matrices  r(l)  to  f  (6)  . 
In  Tables  5  to  6  are  the  magnitudes  and  phases  of  the  coefficients  in 
C  ( 2 )  ^.  Tables  7,  8,  9  and  10  contain  the  magnitudes  of  the  coefficients 
inbi^B(2),  the  difference  C(2)-C(3),  C(4),  and  C(6) ,  respectively.  Some 
notable  observations  are  given  below. 

1.  In  Table  3,  many  gain  estimates  are  less  than  0.8  or  larger  than 
1.2.  In  Table  4,  many  phase  estimates  are  less  than  -10°  or  more 
than  10°.  In  Table  5,  four  coefficients  in  the  main  diagonal  are 
less  than  0.8.  In  Table  6,  three  phase  estimates  in  the  main 
diagonal  deviate  from  zero  by  more  than  10°. 

Remark:  It  appears  that  there  were  difficulties  in  the  initial 

calibration  of  the  antenna  array. 
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2.  In  Table  3,  the  gain  estimates  from  Files  2  and  3  are  essentially 
the  same  and  differ  by  0.034  or  less.  In  Table  4,  with  the  exception 
of  element  #7,  the  phase  estimates  from  these  files  differ  by  1°  or 
less  . 

Remark:  Files  2  and  3  were  collected  with  the  same  source 

elevations  and  ranges.  The  above  difference  in  the  estimates  of 
gains  and  phases  indicates  that  measurement  errors,  including  noise 
and  interference  effects,  were  relatively  small. 

3.  In  Tables  3  and  4,  the  gain  and  phase  estimates  calculated  with  File 
6  differ  noticeably  from  those  calculated  with  Files  2  and  3. 

Remark:  The  result  is  not  unexpected.  It  was  noted  in  Chapter  4 
that  there  are  large  measurement  errors  in  File  6  and  that  these 
errors  cannot  be  suppressed  by  data  preprocessing. 

4.  In  Table  7,  the  magnitudes  of  many  off-diagonal  coefficients  are 
larger  than  0.1. 

Remark:  This  property  indicates  strong  coupling  among  the  array 

elements  or  between  the  elements  and  local  objects,  or  both. 

5.  In  Table  7,  the  magnitudes  for  (p,q)=(l,2)  and  (2,1)  are  0.238  or 
larger.  Those  for  (p,q)=(7,8)  and  (8,7)  are  significantly  smaller 
and  are  0.041  or  less. 

Remark:  The  results  for  (p,q)=(l,2)  and  (2,1)  are  expected.  From 

Table  1,  elements  #1  and  #2  are  separated  by  6.53m,  which  is 
approximately  0.11  wavelength  at  5 . 1  MHz.  With  this  separation,  one 
expects  very  strong  coupling  between  the  two  elements.  Those  for 
(p,q)  =  (7,8)  and  (8,7)  are  unexpected,  because  elements  #7  and  #8  are 
separated  by  6.53m  also.  The  significantly  smaller  values  indicate 
very  strong  interactions  between  the  elements  and  local  objects. 

6.  In  Table  7,  many  off-diagonal  coefficients  in  the  lower  left  corner 
are  larger  than  0.1  while  the  diagonally  opposite  ones  in  the  upper 
right  corner  are  less  than  0.02.  The  coefficients  with  p>9  and  q<;4 
(in  the  lower  left  corner)  are  noticeable  larger  than  those  with  p<;4 
and  qs9  (upper  right  corner) . 
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Remark:  The  true  reason  is  unknown.  These  properties  may  have 

something  to  do  with  the  absence  of  ground  radial  wires  in  elements 
#9  to  #12.  Computer  simulations  show  that  one  can  generate  a  matrix 
bxlB  with  this  property  by  introducing  strong  interactions  between 
these  elements  and  local  objects. 

7.  In  Table  8,  the  magnitude  of  first  coefficient  in  the  main  diagonal 
of  the  difference  C(2)-C(3)  is  zero.  The  magnitudes  of  the  other 
coefficients  are  between  0.012  and  0.025. 

Remarks:  The  first  coefficient  is  zero,  because  C(2)  and  C(3)  both 
satisfy  condition  (5.6).  The  other  coefficients  indicate  that  the 
errors  in  the  main  diagonal  of  Table  5  are  small  and  are  unlikely 
to  exceed  0.03. 

8.  In  Table  8,  five  coefficients  in  the  lower  left  corner,  where  pa9 
and  q<4,  are  larger  than  0.020. 

Remark:  These  values  are  rather  large  compared  with  the  other  off- 
diagonal  coefficients.  They  may  be  related  to  the  absence  of  ground 
radial  wires  for  elements  #9  to  #12. 

9.  The  values  in  the  upper  right  corner  of  Table  10  are  noticeably 
larger  than  those  in  Table  9 . 

Remark:  The  results  are  expected.  The  reasons  are  the  same  as  those 
in  Observation  #3 . 
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Element  Gain  Estimate 

Element 

File=l 

2 

3 

4 

5 

6 

(\J/= 4.5°) 

(ii°) 

(n°) 

(18°) 

(31°) 

(49°) 

1 

1 . 000 

1 . 000 

1 . 000 

1 . 000 

1 . 000 

1 . 000 

2 

1.352 

1.430 

1.408 

1.293 

1.433 

0 . 918 

3 

0 . 870 

0.922 

0 .888 

1.672 

1.218 

4 

1 . 187 

1.261 

1.249 

1.240 

1.691 

1 . 958 

5 

1 . 218 

1.235 

1.225 

1 . 133 

1 . 621 

0 . 810 

6 

1 . 105 

1 .200 

1 . 180 

1 . 116 

1 . 177 

0 .663 

7 

1 . 088 

1.139 

1 . 128 

0 . 978 

1.407 

0 .410 

8 

1 . 127 

1.231  ! 

1.210 

1 . 222 

1.143 

9 

0 . 988 

1.200 

1 . 183 

1.215 

1.611 

0 . 895 

10 

1.079 

1.093 

1.071 

1.064 

11 

1.116 

1.099 

1.134 

1.498 

0 . 861 

12 

1.012 

1.030 

1.013 

0.998 

1.258 

0 . 773 

Table  3:  Estimates  of  element  gains. 
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Element  Phase  Estimate  (deg) 

Element  File=l  I  2  I  3  I  4  I  5  I  6 


(M-5°)  (11°)  (11°)  (18°)  (31°)  (49°) 


9 

8 

10 

11 

16 

8 

-7 

10 

0 

-1 

0 

7 

0 

-11 

11 

-21 

-22 

-21 

-15 

-23 

-34 

12 

-4 

-5 

-5 

1 

-4 

10 

Table  4 .  Estimates  of  element  phases  in  degrees. 


32 


7 . 0 


EVALUATION  OF  CALIBRATION  METHOD 


Matrix  C(2)  ,  the  calibration  matrix  calculated  with  File  2,  was  used 
to  calibrate  the  array  snapshots  in  the  other  files.  The  effect  of 
calibration  on  the  performance  of  the  antenna  array  was  then  studied. 
This  study  included  the  antenna  array  beam  patterns,  the  accuracy  of 
direction  estimates,  and  the  resolution  of  signal  directions. 

It  is  preferable  to  use  C(2)  or  C(3)  to  calibrate  the  array 
snapshots  and  we  have  arbitrarily  chosen  C(2)  .  Matrices  C(2)  and  C(3)  are 
essentially  the  same  and  are  generated  with  Files  2  and  3,  respectively. 
In  Table  2,  these  files  have  the  same  nominal  transmitter  elevations  and 
ranges  (radii) .  Therefore,  we  get  the  best  performance  improvements  the 
calibration  method  can  provide  when  C(2)  is  used  to  calibrate  File  3  or 
C  ( 3 )  is  used  to  calibrate  File  2. 

File  6  was  included  in  the  evaluation.  However,  we  avoid  making 
comments  on  the  positive  results,  if  any,  obtained.  As  noted  at  the 
beginning  of  Chapter  4,  these  results  are  unreliable,  because  of  large 
measurement  errors  associated  with  the  short  transmitter  range  and  the 
high  elevation  angles. 

7.1  Effect  on  Antenna  Beam  Patterns 

The  one-way  ideal  measured  and  calibrated  beam  patterns  of  the 
antenna  array  are  studied  here.  A  one-way  pattern  shows  the  dependence 
of  the  array  output  power  on  signal  source  direction  (6s,\ps)  ,  given  that 
the  source  is  in  the  farfield  zone  of  the  array,  and  the  signal  power  at 
the  array  elements  is  kept  independent  of  the  direction  (6s,\ps).  This 
pattern  depends  on  the  beam  steering  direction,  denoted  by  {6B,\pB)  ,  and 
is  usually  normalized  to  unit  power  at  the  global  maximum.  An  ideal  beam 
pattern  is  a  pattern  produced  with  the  ideal  array  snapshots  given  by 
(5.1)  . 


For  the  sake  of  convenience  in  comparing  beam  patterns,  this  study 
chose  beam  elevation  angles  that  were  equal  to  the  nominal  source 
elevations  in  the  data  files  used  to  produce  the  patterns.  With  these 
choices,  there  is  some  assurance  that  the  main  peaks  of  the  ideal, 
measured  and  calibrated  beam  patterns  would  be  close  to  the  beam  steering 
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directions . 


The  calculation  of  beam  patterns  must  take  into  consideration  that 
the  signal  power  in  the  array  snapshots  is  not  fixed  in  every  data  file. 
There  are  two  reasons.  Firstly,  as  one  can  observe  from  Figs.  3  to  8,  the 
aircraft  flight  paths  were  not  perfect  circles  about  the  vertical  W-axis 
of  the  UVW- coordinate  system.  Deviations  from  the  circular  path  produced 
variations  of  signal  power  at  the  array  elements.  Secondly,  because  of 
air  turbulence  behind  an  aircraft,  the  transmitter  orientation  was  not 
fixed  relative  to  the  aircraft.  Changes  in  transmitter  orientations  also 
produced  variations  in  signal  powers  at  the  array  elements. 

Here,  we  attempt  to  suppress  the  dependence  of  signal  power  on 
source  direction  by  normalizing  separately  the  measured  and  calibrated 
array  snapshots,  which  are  complex  vectors  with  N  components,  to  vectors 
with  unit  magnitudes  before  using  them  in  the  calculations. 

The  three-steps  procedure  described  below  was  used  to  calculate  the 
beam  patterns.  In  Step  2,  M  is  the  number  of  array  snapshots  in  a 
preprocessed  file,  the  set  {  (01(  </y)  ,  ( d2 ,  \p2)  ,  •  •  •  ,  (0M,  \pn)  }  denotes  the  source 
directions,  {x1 ,  x2 ,  •  •  •  ,  xM}  denotes  the  preprocessed  snapshots,  and 
{x^,x2,  •  •  •  ,x^}  denotes  the  corresponding  calibrated  snapshots.  These 
calibrated  snapshots  were  calculated  with  (5.14)  and  the  matrix  C(2). 

Step  1 . 

Specify  the  beam  steering  direction  {dB,\pB}  . 

Denote  the  corresponding  array  steering  vector  by  a{6B,\pB)  . 

Step  2 . 

For  m=l , 2 , • • • , M, 

(a)  Calculate  two  unit -magnitude  vectors  x^  and  x^  as 


K  =  /Kl  -  n.i) 

K  =  K  /Kl  '  {1.2) 

(b)  Calculate  the  beam  amplitudes  in  the  ideal,  measured  and 
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calibrated  beam  patterns  as 


Bmeas(0B^B;0m^J  =  |aH(0B,^B)im|2  , 

Bcai^s'Vv^’U  =  laH(0B'^B)^l2  ' 


Step  3 

(a)  Find  the  largest  amplitude  in  the  set  {Bideal  (  0B,  ^B;  0m,  V'J  : 
m=l ,  2 ,  •  •  •  ,  M}  ,  which  is  calculated  with  (7.3)  .  Use  it  to  divide 
every  amplitude  in  this  set. 

(b)  Repeat  with  the  other  two  sets  calculated  with  (7.4)  and 
(7.5)  . 

Figs.  13  and  14  to  17  contain  some  results  that  are  typical  of  the 
beam  patterns  obtained  with  Files  1  and  3  to  6,  respectively.  They  are 
obtained  with  the  beam  elevation  \pB  equated  to  the  nominal  source 
elevations  in  Table  2  and  the  beam  bearing  arbitrarily  assigned  as  0B=45°. 
Each  figure  is  divided  into  two  halves:  the  upper  half  shows  the  patterns 
for  all  source  bearings,  and  the  lower  half  is  a  magnification  of  the 
section  with  0°<6S<90°  .  Vertical  lines  have  been  drawn  through  the  source 
bearing  6S=6B to  facilitate  the  study  of  the  patterns  in  this  region.  The 
figures  show  the  dependence  of  beam  amplitudes  on  source  bearing  6S  for 
0b=45°  and  \pB=\ps=4 . 5°  ,  11°  ,  18°  ,  31°  ,  49°  .  Some  notable  observations  are  given 
below : 

1.  The  main  peaks  of  the  ideal  beam  patterns,  drawn  in  solid  lines,  in 
the  upper  halves  of  Figs.  13  to  17  are  centred  at  0S=45°,  the 
bearing  of  the  beam  steering  direction.  Up  to  eleven  sidelobes  are 
present.  Some  gradually  merge  together  as  the  beam  elevation 
increases.  The  peak  values  of  many  sidelobes  are  rather  high  and  are 
larger  than  -10  dB . 

Remark:  The  above  properties  of  the  ideal  beam  patterns  are 

determined  by  the  design  of  the  antenna  array.  They  have  nothing  to 
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POWER  (DB)  POWER  (DB) 


SOURCE  BEARING  C  DEG ) 


SOURCE  BEARING  ( DEG ) 


Fig.  13  The  ideal  (solid  line),  measured  (dotted)  and  calibrated 
(broken)  array  beam  patterns  for  beam  steering  direction  (0,f)= 
(45°, 4. 5°)  and  signal  source  elevation  f e=4 . 5 .  The  lower  figure  is  a 
magnification  of  the  section  with  source  bearing  0°  to  90°  and  the 
vertical  lines  are  drawn  through  the  beam  bearing  0=45°  .  File  1  has  been 
used  to  produce  the  figure. 
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POWER  (DB)  POWER  (DB) 


SOURCE  BEARING  (DEG) 


SOURCE  BEARING  (DEG) 


Fig.  14  The  ideal  (solid  line),  measured  (dotted)  and  calibrated 
(broken)  array  beam  patterns  for  beam  steering  direction  (0,f)= 
(45°,  11°).  File  3,  where  fs=ll°,  has  been  used. 
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POWER  (DB)  POWER  (DB) 


SOURCE  BEARING  (DEG) 


SOURCE  BEARING  ( DEG ) 


Fig.  15  The  ideal  (solid  line),  measured  (dotted)  and  calibrated 
(broken)  array  beam  patterns  for  beam  steering  direction  (8,f)= 

(45°, 18°).  File  4,  fs=18°,  has  been  used. 
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POWER  (OB)  POWER  (OB) 


Fig.  16  The  ideal  (solid  line),  measured  (dotted)  and  calibrated 
(broken)  array  beam  patterns  for  beam  steering  direction  (0,f)= 
(45° ,31°).  File  5,  fs=31°,  has  been  used. 
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POWER  (DB) 


SOURCE  BEARING  (DEG) 


SOURCE  BEARING  (DEG) 


Fig.  17  The  ideal  (solid  line),  measured  (dotted)  and  calibrated 
(broken)  array  beam  patterns  for  beam  steering  direction  (0,t)= 
(45°, 49°).  File  6,  fs=49°,  has  been  used. 
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do  with  the  calibration  experiment. 


2.  The  measured  beam  patterns,  in  dotted  lines,  look  similar  to  the 
ideal  ones.  There  are  biases  in  the  positions  of  the  main  peaks, 
because  they  are  not  centred  at  0S=45°.  These  biases  are  easier  to 
observe  when  the  lower  halves  of  the  figures  are  examined. 

Further  studies,  to  be  presented  later,  show  that  the  biases  depend 
on  source  directions. 

Remark:  The  similarity  indicates  that  the  array  snapshots  can  be 

used  to  estimate  signal  directions.  The  biases  indicate  that  these 
estimates  require  corrections  that  depend  on  source  directions. 

3.  The  calibrated  beam  patterns,  in  broken  lines,  also  look  similar  to 
the  ideal  ones.  In  Figs.  13  to  15,  where  \pB=\ps&18  ° ,  the  amplitudes 
in  the  three  strongest  lobes  are  closer  to  the  ideal  patterns  than 
the  measured  ones.  In  Figs.  16  and  17,  where  '/'B=31  and  49°, 
respectively,  the  identifiable  effect  of  calibration  is  only  a 
reduction  of  the  biases  in  the  main  lobe  positions. 

Remark:  The  calibration  matrix  calculated  with  File  2,  where  the 

source  elevation  is  11°,  is  effective  for  calibrating  Files  1,  3  and 

4.  It  is  ineffective  for  Files  5  and  6. 

4.  In  a  separate  study,  not  present  here,  the  calibration  matrix 
calculated  with  File  6  was  used  to  calibrate  File  6 .  The 
identifiable  effect  of  this  calibration  was  also  a  reduction  of  the 
bias  in  the  main  lobe  position  only. 

Remark:  The  calibration  method  cannot  be  used  to  calibrate  File  6, 
because  it  has  large  measurement  errors. 

7.2  Effect  on  Direction  Estimates 

This  study  is  restricted  to  the  special  case  where  there  is  only  one 
signal  in  each  array  snapshot  and  that  this  signal  is  the  direct  signal 
from  the  transmitter.  It  uses  the  parametric  target  model  fitting  ( PTMF) 
method  [13-15] ,  which  is  an  HRDF  method,  to  estimate  the  source  bearing 
and  elevation  in  the  snapshots. 
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In  the  presentation  of  results,  the  error  in  a  bearing  estimates  8 
is  defined  as 


where  8  is  the  GPS  bearing  calculated  with  (3.2)  and  (3.3). 

The  errors  in  bearing  estimates  are  plotted  in  point  mode  in  Figs. 
18  to  22.  Those  of  the  elevation  estimates  are  in  Figs.  23  to  27, 
respectively.  Some  notable  observations  are  given  below. 

1.  The  estimates  seem  to  fall  on  a  smooth  curve,  although  they  are 
plotted  in  point  modes . 

Remark:  It  appears  that  the  SINRs  in  the  preprocessed  array 

snapshots  are  high  and  that  the  interfering  environment  did  not 
change  rapidly  with  source  bearings. 

2.  Before  calibration,  the  errors  in  bearing  estimates  depend  strongly 
on  the  source  bearings.  Generally,  they  are  positive  in  the  60° 
region  and  are  negative  in  the  240°  region.  The  calibration 
suppresses  the  dependence  of  errors  on  source  bearings . 

3 .  The  calibration  reduces  the  average  magnitude  of  errors  in  bearing 
estimates.  Many  reductions  are  more  than  one  degree  in  Figs.  18  to 
20,  where  the  source  elevations  are  18°  or  less. 

4.  There  are  strings  of  zero  elevation  estimates  in  Fig.  23,  where  the 
source  elevation  is  4.5°. 

Remark:  The  elevation  estimates  obtained  from  File  1  are  not 

expected  to  be  accurate.  In  this  file,  the  antenna  array  aperture 
in  the  vertical  direction  is  less  than  0.2  wavelength  for  all  source 
bearings.  This  aperture  is  too  small  to  produce  reliable  elevation 
estimates . 

Remark:  A  separate  study  showed  that,  when  an  elevation  estimate 
was  zero,  the  PTMF  method  could  usually  produce  a  better  fit  to  the 
data  if  it  was  allowed  to  use  a  value  of  cos(i })  that  was  larger  than 
unity.  This  cos(i })  corresponded  to  a  non-physical  value  of  i p . 
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Bearing  (deg) 


Figure  18.  Errors  in  bearing  estimates  obtained  with  File  1.  The  upper 
figure  is  calculated  with  the  uncalibrated  data  and  the  lower  figure  is 
calculated  with  the  calibrated  data. 
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Bearing  Error  e(0)  (deg)  Bearing  Error  e(0)  (deg) 


5 


Figure  20.  Errors  in  bearing  estimates  obtained  with  File  4. 


Figure  21.  Errors  in  bearing  estimates  obtained  with  File  5. 
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Bearing  Error  e(0)  (deg)  Bearing  Error  e(0)  (deg) 
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Elevation  (deg)  Elevation  (deg) 


Figure  23. 


Elevation  estimates  obtained  with  File  1 


Elevation  (deg) 


30 


Calibrated 


Elevation  (deg)  Elevation  (deg) 
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Uncalibrated 


30  h 


Bearing  (deg) 


Figure  25.  Elevation  estimates  obtained  with  File  4. 


56 


Elevation  (deg)  Elevation  (deg) 


Bearing  (deg) 


Figure  26.  Elevation  estimates  obtained  with  File  5. 
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Elevation  (deg)  Elevation  (deg) 


Figure  27 .  Elevation  estimates  obtained  with  File  6 . 
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5.  The  upper  halves  of  Figs.  2  4  to  26  show  that  the  elevation  estimates 
from  Files  3  to  5,  respectively,  depend  on  the  source  bearings.  In 
particular,  they  are  below  the  nominal  values  in  the  neighbourhood 
of  bearing  0S=13  5°  and  above  the  nominal  values  in  the  neighbourhood 
of  <9s  =  315°.  This  dependence  is  not  visible  after  calibration. 

6.  There  are  biases  in  the  elevation  estimates  in  Figs.  24  to  27.  They 
depend  on  the  files  and  are  more  noticeable  in  the  lower  halves  of 
the  figures. 

Remark:  There  are  many  possible  reasons  for  these  biases.  They 

include  flight  path  distortions,  the  deviation  of  transmitter 
altitudes  from  nominal  altitudes,  and  the  dependence  of  the  A'  in 
(3.2)  on  data  files. 

7.  In  both  the  upper  and  lower  halves  of  Fig.  27,  there  are  sudden 
drops  of  about  10°  in  elevation  estimates  at  0s=172°. 

Remark:  After  preprocessing,  the  initial  and  final  transmitter 

bearings  in  File  6  are  172°  and  171.9°,  respectively.  From  Fig.  8, 
the  transmitter  ranges  are  different  for  these  two  bearings.  It  is 
highly  probable  that  the  actual  transmitter  altitudes  are  also 
different . 

7.3  Effect  on  Resolution  of  Signal  Directions 

This  study  compares  the  ability  of  the  PTMF  method  to  resolve  the 
directions  of  two  signal  sources  before  and  after  calibration.  It 
involves  the  construction  of  array  snapshots  with  two  signal  sources 
present . 

The  study  uses  the  five-step  procedure  described  below.  In  Step  1, 
a  section  of  2K  preprocessed  array  snapshots,  denoted  by 
{ x ” , x " , x " ,  *  * ’ , x 2 K }  ,  is  extracted  from  a  specified  file.  This  section  is 
60°  wide  in  signal  source  bearing  and  the  source  bearing  at  the  centre 
of  the  section  is  close  to  30°.  In  Step  2,  this  section  is  separated  into 
two  groups  given  by  {x'^x'^x'J,  •  *  •  and  {xV,K,  x'^K_2  ,  x’^K_4 ,  *  •  •  ,x’’}  .  Because 

the  ordering  of  the  snapshots  in  the  second  group  has  been  reversed,  the 
source  bearing  increases  from  0°  to  60°  in  one  group  and  decreases  from 
60°  to  0°  in  the  other.  These  two  groups  are  then  combined  to  produce  a 
new  set,  denoted  by  {Yi  #  Y2  >  Yi '  *  ’  *  '  Yk)  '  ^-n  which  there  are  two  signal 
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sources  crossing  each  other  near  the  30°  bearing.  In  Step  4,  (7.8)  is 
obtained  from  (5.14)  by  using  C(2)  instead  of  C,  yk  instead  of  x,  and  yk 
instead  of  xc.  The  details  are  given  below: 

Step  1 . 


Extracts  a  section  of  2K  snapshots  with  source  bearings  between  0° 
and  60°  from  a  specified  file.  Identify  the  snapshots  by 
{X'l'X'2,X'3,-.-,X'2K}. 

Step  2 . 


Construct  a  set  of  K  new  snapshots  {y1,y2,y3,  •  ■  •  ,yK}  as 


Y>.  =  X2k-i  +  x: 


2K+2 -  2k  ' 


k=l , 2 , 3 , .  .  .  ,K  . 


(7.7) 


Step  3  . 


Use  the  PTMF  method  to  estimate  two  signal  directions  from  each 
snapshot  in  {y1,y2,y3,  •  •  -,yK}  .  Denote  the  estimates  obtained  from  yk 
bY  and  {Sk2,ik2}. 

Step  4  . 


Construct  a  set  of  K  calibrated  snapshots  {y3,y2,y3 


/Yk) 


as 


Yk  =  C(2)yk 


k=l , 2 , 


i  K 


(7.8) 


Step  5  . 


Use  the  PTMF  method  to  estimate  two  signal  directions  from  each 
snapshot  in  {y3 ,  y\ ,  y3 ,  •  •  •  ,  y^}  .  Denote  them  by  {Ki'Ki)  and  {^2-^2}- 

The  value  of  K  is  300  for  Files  1  and  3  to  5,  and  170  for  File  6. 
It  is  smaller  for  File  6,  because  this  file  is  smaller  and  the  array 
snapshots  are  spaced  farther  apart  in  bearing. 

Figs.  28  to  32  contain  the  bearing  estimates  0kl,  6k2,  6ckl,  and  0k2 
obtained.  These  estimates  are  plotted  in  point  mode.  Each  figure  has  two 
solid  lines.  They  are  nominal  source  tracks  calculated  with  (3.2)  and 
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Bearing  Estimates  (deg)  Bearing  Estimates  (deg) 


Snapshot  # 


Figure  28.  Bearing  estimates  (points)  calculated  with  File  1.  The  solid 
lines  represent  the  nominal  source  tracks. 
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Bearing  Estimates  (Deg)  Bearing  Estimates  (Deg) 


Snapshot  # 


Snapshot  # 

Figure  29.  Bearing  estimates  (points)  calculated  with  File  3. 
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Bearing  Estimates  (Deg)  Bearing  Estimates  (Deg) 


Snapshot  # 

Figure  30.  Bearing  estimates  (points)  calculated  with  File  4. 
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Bearing  Estimates  (Deg)  Bearing  Estimates  (Deg) 


uncalibrated 


Snapshot  # 


Figure  31.  Bearing  estimates  (points)  calculated  with  File  5 


Bearing  Estimates  (Deg)  Bearing  Estimates  (Deg) 


Snapshot  # 


Figure  32.  Bearing  estimates  (points)  calculated  with  File  6.  The  true 
signal  tracks  are  given  by  the  solid  lines. 
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(3.3).  One  is  the  source  track  in  {x» ,  x" ,  x'^,  •  •  •  ,  x'^ }  .  The  other  is  the 

source  track  in  {x2K'x2k-2'x2k-4'  ‘  ‘  '  »x"}  •  There  are  several  notable 

observations : 

1.  In  each  figure,  the  bearing  estimates  seem  to  cluster  around  two 
apparent  tracks  that  intersect  near  the  30°  bearing. 

2 .  The  apparent  tracks  deviate  from  the  nominal  tracks  given  by  the 
solid  lines.  In  particular,  all  are  above  the  nominal  tracks  before 
calibration. 

Remark:  The  deviations  are  related  to  the  errors  in  bearing 

estimates  in  Figs.  18  to  22.  Before  calibration,  the  bearing 
estimates  are  larger  than  the  nominal  bearings  in  the  bearing  range 
0°  to  60°.  This  agrees  with  the  observation  that  the  apparent  tracks 
are  above  the  nominal  tracks. 

3.  The  scattering  of  the  points  about  the  apparent  tracks  is  generally 
larger  if  the  angular  separation  of  the  true  signal  bearings  is 
smaller. 

Remark:  This  property  is  expected,  because  it  is  more  difficult  to 
resolve  two  signal  directions  that  are  closer  together  in 
directions.  More  difficulty  in  resolution  leads  to  larger  errors  in 
the  estimates. 

4.  The  calibration  reduces  the  scattering  of  the  bearing  estimates 
about  the  apparent  tracks. 
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8.0  CONCLUSIONS 


The  authors  have  carried  out  a  second  analysis  of  the  5.1  MHz 
calibration  data  measured  with  a  high-frequency  (HF)  antenna  array  at 
Leitrim.  This  analysis  used  the  second  antenna  array  calibration  method 
[1]  developed  by  Eric  Hung  at  DREO.  It  was  part  of  a  task  to  develop  a 
calibration  method  for  the  HF  antenna  arrays  in  the  Canadian  Forces 
Supplementary  Radio  Systems  (CFSRS)  modernization  project.  These  arrays 
are  used  for  detecting  long-range  HF  signals  and  estimating  their 
directions . 

The  following  observations  are  based  on  the  analysis  of  data  files 
generated  with  signal  source  elevations  31°  or  lower: 

1.  The  calibration  method  reduces  the  deviation  of  the  antenna  beam 
patterns  from  the  ideal  patterns; 

2 .  The  method  increases  the  accuracy  of  direction  estimates  and 
suppresses  the  dependence  of  errors  in  bearing  estimates  on  source 
bearings;  and 

3.  The  method  increases  the  ability  of  an  HRDF  method  to  resolve  signal 
directions . 

4.  The  method  is  superior  to  the  first  antenna  array  calibration 
method.  In  particular,  the  direction  estimates  are  more  accurate  and 
the  ability  of  an  HRDF  method  to  resolve  signal  directions  is 
better . 

It  is  not  possible  to  remark  on  the  results  obtained  with 
calibration  data  generated  with  signal  source  elevations  higher  than  31°  . 
The  results  are  very  encouraging.  However,  these  data  files  have  large 
errors  that  cannot  be  suppressed. 
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